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The relationship between the pionium lifetime and the tttt scattering lengths is established, including the sizable 
electromagnetic corrections. The bound state formalism that is used is that of constraint theory which provides 
a covariant three-dimensional reduction of the Bethe-Salpeter equation. The framework of generalized chiral 
perturbation theory allows then an analysis of the lifetime value as a function of the tttt scattering lengths, the 
latter being dependent on the quark condensate value. 



The possible measurement of the pionium 
(tt+tt" atom) lifetime with a 10% precision in the 
DIRAC experiment at CERN is expected to al- 
low a determination of the combination (oq — Cq) 
of the TTTT scattering lengths with 5% accuracy; 
here, Oq is the strong interaction (dimensionless) 
S-wave scattering length in the isospin / chan- 
nel. The strong interaction scattering lengths Oq 
and Oq have been evaluated in the literature in the 
framework of chiral perturbation theory [xPT) to 
two-loop order of the chiral effective lagrangian 
|2|-Q. Therefore, the pionium lifetime measure- 
ment provides a high precision experimental test 
of chiral perturbation theory predictions. 

The nonrelativistic formula of the pionium life- 
time in lowest order of electromagnetic interac- 
tions was first evaluated by Deser et al. It 
reads: 
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where Am^ = — m^o and tp^ (0) is the wave 

function of the pionium at the origin (in a;-space) . 

A precise comparison of the theoretical values 
of the strong interaction scattering lengths with 
experimental data necessitates, however, an eval- 
uation of the corrections of order 0{a) (a being 
the fine structure constant) to the above formula. 
Such an evaluation was recently done by several 
authors. In the frameworks of quantum field the- 



ory and xPT, three different methods of evalu- 
ation have led to the same estimate, of the or- 
der of 6%, of these corrections |^-^. The first 
method uses a three-dimensionally reduced form 
of the Bethe-Salpeter equation (constraint theory 
approach) and deals with an off-mass shell formal- 
ism 1^,0. The second method uses the Bethe- 
Salpeter equation with the Coulomb gauge 
The third one uses the approach of nonrelativis- 
tic effective theory [||. 

The pionium lifetime, with the sizable 0{a) 
corrections included in, can be represented as: 
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where 7?.eA4oo,H is the real part of the on- 
mass shell scattering amplitude of the process 
TT+TT^ 7r°7r°, calculated at threshold, in the 
presence of electromagnetic interactions and from 
which singularities of the infra-red photons have 
been appropriately subtracted |l^; the factor 7 
represents contributions at second-order of per- 
turbation theory with respect to the nonrela- 
tivistic zeroth-order Coulomb hamiltonian of the 
bound state formalism. The explicit expressions 



of 7?.eA^oo,H s-nd of 7 may differ from one ap- 
proach to the other, but their total contribution 
should be the same. 

Our aim is to extend the previous analysis to 
the case of generalized llU'S ■ '^^^ latter is 

based on the observation that the fundamental 
order parameter of spontaneous chiral symmetry 
breaking is -F^ , the decay coupling constant of the 
pion, which is related to the two-point function of 
left- and right-handed currents in the chiral limit. 
The other order parameters, such as the quark 
condensate in the chiral limit, < OjgqlO >o, have 
values depending on the details of the mechanism 
of chiral symmetry breaking and require indepen- 
dent experimental tests. Standard xPT is based 
on the assumption that the value of the Gell- 
Mann-Oakes-Renner (GOR) parameter [Q, de- 
fined as 
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where 2m = to„ -|- and is the pion physical 
mass, is close to one. Stated differently, the quark 
condensate parameter 
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where F is in the chiral SU{2) x SU{2) limit, 
is of the order of the hadronic mass scale Ajf ~ 1 
GeV. This assumption fixes the way standard 
xPT is expanded: the quark condensate parame- 
ter Bq is assigned dimension zero in the infra-red 
external momenta of the Goldstone bosons, while 
the quark masses are assigned dimension two |^] . 

Generalized relaxes the previous assump- 
tion and treats the order of magnitude of the 
quark condensate parameter Bq as an a priori 
unknown quantity (awaiting a precise experimen- 
tal information about it) leaving to it the possibil- 
ity of reaching small or vanishing values. To this 
aim, Bq is assigned dimension one in the infra- 
red momenta of the external Goldstone bosons 
and accordingly quark masses are also assigned 
dimension one. Due to this rule, at each order 
of the perturbative expansion, generalized x^T 
contains more terms than standard xPT . For in- 
stance, the pion mass formula becomes, at leading 
order, = 2rhBo + Am^Ao, where the constant 



Aq is expressible in terms of two-point functions 
of scalar and pseudoscalar quark densities. In the 
standard xPT case, this term is relegated to the 
next-to-leading order. 

At the tree level of the chiral effective la- 
grangian, the tttt scattering amplitude A(s\t,u) 
has the expression A{s\t,u) = (s — 2mBo)/F'^, 
which displays explicit dependence on the quark 
condensate parameter. It is useful to introduce 
two parameters, a and /3, that allow one to ex- 
press the amplitude A in terms of the physical 
constants Ft^ and rriTr OjF 
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At leading order, F^^ ~ F and f3 — 1. The 
GOR parameter (H) is then related to the pa- 
rameter a by the relation: xgor = (4 — a)/3. 
For a = 1, one recovers the standard xPT case, 
where xgor = 1- When a increases, xgor de- 
creases, corresponding to decreasing values of the 
quark condensate. For a = 4, one reaches the ex- 
treme case of generalized xPT, where xgor and 
the quark condensate vanish. At higher orders, 
the relationship between a and xgor becomes 
more complicated, involving the parameter f3 and 
other low energy constants, but the main quali- 
tative feature found above remains: values of a 
close to 1 correspond to the case of standard xPT, 
while large values of a, above 1.5-2, say, cover the 
complementary domain of generalized xPP ■ 

At two-loop order, the tttt scattering amplitude 
is described by six parameters, a, (3, Ai, A2, A3 
and A4. The five parameters other than a are 
weakly dependent on the quark condensate and 
the values of the A's are fixed by a detailed anal- 
ysis of sum rules using available data on tttt scat- 
tering at medium energies. The parameter (3 is 
essentially sensitive to the deviation of from 
F and remains in general close to 1. One ends 
up, in generalized xPT ^ with a complete deter- 
mination of the threshold parameters of the tttt 
scattering amplitude in terms of the possible nu- 
merical values of the parameter a together with 
a best value of the other parameters. The expres- 
sions of the scattering lengths to two-loop order, 
as well as the values of the parameters A, can be 
found in Rcf. [||. 



The evaluation of the corrections contained in 
the decay width formula (||) can be done in much 
the same way as in the standard case, with 
the difference that whenever scattering lengths 
or the quark condensate parameter appear, these 
should be expressed through the generalized xPT 
formulas or parametrizations. We have followed 
the same method of approach as in Refs. 
and skip here the details of the calculations. We 
only note that we made the approximation of tak- 
ing the same values for the electromagnetic low 
energy constants ki in generalized and standard 
XPT Actually, the effective lagrangian 

with electromagnetism is not yet available in gen- 
eralized xPT, where the number of such con- 
stants should be larger. However, the contri- 
butions of the constants ki not being dominant 
in the total decay width correction, such an ap- 
proximation seems to be justified. Furthermore, 
a 100% uncertainty has been assigned to their 
global contributions in the various types of cor- 
rection, their values being taken from Ref. [p^ . 
Also, weakly dependent terms on the parameter 
a have been numerically taken the same as in the 
standard case. Numerical estimates have been 
done with the following inputs: F-^ = 92.4 MeV, 
= 139.57 MeV, m^o = 134.97 MeV. 

We have chosen to analyze the total correction 
in terms of the scattering lengths, which are the 
physical quantities of interest, rather than of the 
intermediate parameters a and (5. Since the pio- 
nium lifetime concerns at leading order the com- 
bination (aq — Cq) of the scattering lengths, it 
is natural to promote this quantity as the main 
variable of the problem. As a second variable, 
we choose the P-wave scattering length, which 
is most sensitive to the parameter (5 and almost 
insensitive to a. One then expresses the total cor- 



rection in terms of (oq - 



and the A's. (It is 



sufficient for the correction term to use the one- 
loop expressions of the scattering lengths, from 
which A3 and A4 are absent.) The experimental 
value of a\m\ is 0.038 ± 0.002 ll6|. One can then 
study the dependence of the decay width (or of 



the lifetime) upon the combination (oq 



of 



the scattering lengths for fixed values of the above 
parameters. 

The previous analysis can be repeated with var- 



ious sets of inputs for a\ and for the A's within 
the intervals of their possible values. Generally, 
for a given (og — Oq), the lifetime t varies little for 
such changes of the inputs. The variations of the 
A's induce at most a variation of 0.8% on r, while 
the variations of a\ induce at most a variation of 
1.1%. Assuming these variations as being uncor- 
related, one may consider that 2% is an upper 
bound for the possible variations of r. Consider- 
ing the latter as un uncertainty and adding it to 
the 2% uncertainty obtained in the course of cal- 
culation of At/tq, one obtains a total uncertainty 
of 4% around the value of r calculated with the 
central values of a\ and of the A's. 

The above analyzes are graphically summa- 
rized in Fig. ^ The full line represents the 
lifetime r as a function of the the combination 
(oq — Oq), corresponding to the central values of 
a} and of the A's. The band around it corresponds 
to the estimated 4% total uncertainty. As already 
mentioned, the parameter is almost insensitive 
to the variations of a; for instance, on the central 
line, it varies between 1.10 and 1.11. For a given 
experimental value of the lifetime, with a possi- 
ble uncertainty, one can deduce from Fig. Q the 
corresponding value of the combination (oq — Og) 
with the related uncertainty. 

The interpretation of the experimental value 
of the pionium lifetime will depend on its or- 
der of magnitude. We assume in the follow- 
ing discussion that, as the DIRAC experimen- 
tal project foresees, the experimental uncertainty 
on the lifetime is of 10%. For larger uncertain- 
ties, the discussion should accordingly be mod- 
ified. Three different possibilities may be con- 
sidered. First, the central value of the lifetime 
is close to 3 X 10~^^ s, lying above 2.9 x 10~^^ 
s, say. Then standard is firmly confirmed, 

since its predictions of (og — Og) lie between 0.250 
and 0.258 Second, the central value of the 
hfetimc hes below 2.4 x 10^^^ s. Then it is the 
scheme of generalized xPT which is confirmed, 
since the corresponding central values of a would 
lie above 2, which would mean that the quark 
condensate is not as large as assumed in the stan- 
dard scheme. The third possibility is the most 
difficult to interpret. If the central value of the 
lifetime hes in the interval 2.4 — 2.9 x 10^^^ s. 
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Figure 1. The pionium lifetime, in units of 10^^^ 
s, as a function of the combination {a^ — Oq) of 
the 5*- wave scattering lengths. The full line corre- 
sponds to the central value ajmj = 0.038 of the 
P-wave scattering length. The band delineated 
by the dotted lines takes account of the uncer- 
tainties, coming from theoretical evaluations, low 
energy constants and a\. 



then, because of the uncertainties, ambiguities in 
the interpretation may arise. Let us consider as 
an illustrative possibility the hypothetical result 
r = (2.60 ± 0.26) x lO'^^ s. Without the in- 
clusion of the uncertainties discussed above, this 
would imply (a[] - Oq) = 0.278 =F 0.015 and a = 
1.68 =F 0.42. Taking into account the 4% uncer- 
tainties represented by the band, one would have: 
(aO - al) = 0.278 T 0.022 and a = 1.68 T 0.61. 
Clearly, the upper bound of r is in favor of stan- 
dard xPT, while its lower bound in favor of gen- 
eralized X-PT. Such a situation would not dis- 
entangle the two schemes from each other. The 
uncertainties might be reduced in this case by a 
more refined analysis, taking into account the cor- 
relations between variations of a\ and of the A's. 
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